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Abstract: Non-Fellerian processes show phenomena that are unseen in standard inter¬ 
acting particle systems. We consider freezing transitions in one-dimensional non-Fellerian 
processes which are built from the abelian sandpile additions to which in one case, spin 
flips are added, and in another case, the so called anti-sandpile subtractions. In the first 
case and as a function of the sandpile addition rate, there is a sharp transition from a 
non-trivial invariant measure to the invariant measure of the sandpile process. For the 
combination sandpile plus anti-sandpile, there is a sharp transition from one frozen state 
to the other anti-state ^. 


1 Introduction 

Much of the motivation in the study of interacting particle systems has come from the 
search for new phenomena. The abelian sandpile model has been widely studied in 
the context of so called self-organized criticality. From the point of view of probability 
theory, it is the best known case of a spatially extended non-Fellerian stochastic dynam¬ 
ics, HH. It has challenged our basic understanding of the construction of interacting 
processes in inhnite volume, even in one dimension. 

In one dimension, the stationary measure of the standard abelian sandpile model is 
trivial in the thermodynamic limit. However the dynamics of relaxation to this measure 
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is non-trivial. In |T21 we have constructed the dynamics for the one-dimensional sandpile 
model on the inhnite lattice Z. The result is a monotone non-Fellerian process which 
converges in hnite time to its unique stationary state, which is concentrating on the 
maximal conhguration. As soon as one changes to other lattices, such as decorated 
one-dimensional lattices, the triviality of the limiting stationary measure disappears, 
and especially in one dimension existence of thermodynamic limits is not guaranteed 
due to the presence of inhnite avalanches [H]. 

In the present paper, we combine the one-dimensional sandpile model with a spin- 
hip dynamics (pure spin hip as well as Glauber type or more general spin hip processes 
with positive rates). Indeed, in one dimension, the standard sandpile model has only 
two possible heights per site, and spin hip just means changing the height from one to 
the other. In the language of sandpiles, adding a pure spin hip is the simplest example 
of combining two diherent toppling mechanisms, the spin-hip part corresponding to a 
purely dissipative (diagonal) toppling matrix. More precisely, in Section 2, our dynamics 
has a formal generator 

Lfiv) = « - fiv)] + “ /(^)] 

X 

where denote the abelian sandpile addition operators, and 9x the hip operator, on 
conhgurations rj G {1,2}^. In words that means that at rate a we add and stabilize 
according to the abelian sandpile rule, and at rate c{x,ri), we just hip the value of 
the height where we added. In that way we have a parameter a that describes the 
relative weight of the sand additions versus the spin hips. The resulting “sand-hip” 
dynamics shows a freezing transition as a function of that a. In the simplest case where 
c{x,r]) = 1 (adding pure spin hip), our main result says that for a > 1 there is a 
hnite time after which the system reaches the maximal conhguration (i.e., the sandpile 
part “wins”), whereas the unique stationary measure is non-trivial and mixing under 
spatial translations for a < 1. That is a strong manifestation of the non-locality of the 
dynamics. Indeed, for Fellerian processes such phenomenon cannot occur. 

As a second example of such “competition of diherent additions”, in Section 3 we 
consider a combination of a sandpile and an anti-sandpile process. This dynamics is 
inspired by jHl. The anti-sandpile part of the dynamics consists of removing grains 
and stabilizing by reverse topplings. The inhnite volume limit of the anti-sandpile 
stationary measure is a Dirac measure concentrating on the minimal conhguration. 
Our main result here is that unless the rates of addition and subtraction are equal, the 
limiting stationary measure is a Dirac measure corresponding to the dominant rate. We 
thus have a sharp transition between two diherent frozen states. 

These nonequilibrium phase transitions have an interest of their own but they also go 
some way in adding extra and physically relevant interactions to the standard abelian 
sandpile. We have in mind the sticky sandpiles of [S], for which our dynamics is a 
subclass, and for which various transitions have been numerically checked. 

As a hnal note, it is interesting to make the analogy with non-Gibbsian measures. 
In some sense, they are the “equilibrium analogue” of non-Fellerian interacting particle 
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systems. In ^3] an example of a freezing transition was obtained, strongly connected 
with the absence of continuity of the local conditional probabilities. 

2 Adding spin flips to the sandpile process 

The state space of our process is = {1, 2}^. For a conhguration rj E fl, rj^x) G {1, 2} 
is usually interpreted as the height or the number of grains at site x. That language 
can be continued even when combining the sandpile automaton with other dynamics 
but we prefer to use the words “active” for r]{x) = 2 and “inactive” for r]{x) = 1. 

The dynamics will change the conhguration according to a combination of the standard 
sandpile model and a spin hip dynamics. We start with the simplest form of spin hip, 
changing “active” into “inactive” and vice versa at rate 1: the spin hip 9^ is thus dehned 
as 


{r]{x) + 1) mod 2, ii y = x 
viy)^ iiy 7^ X 


OxViy) 


( 2 . 1 ) 


Only in Section 12.41 will we generalize the spin hip part of the dynamics. 

For the sandpile dynamics, we can rely on our previous work in ^21 where we have 
studied the inhnite volume limit of the one-dimensional sandpile process. We will 
therefore not bother to redo the limiting procedures but below we immediately give 
the result, the form of the inhnite volume addition operators a^- The informal verbal 
prescription of the sandpile dynamics goes as follows: if a site is inactive, it becomes 
active at rate a. If the site x is already active, one looks left and right of x at the 
closest sites x~ and x'^ which are inactive. Again at rate a these two become active 
and the mirror image of x with respect to the middle of [a;“, becomes inactive. That 
corresponds, in the inhnite volume limit, to the result (in hnite volume) of adding and 
stabilizing through a sequence of topplings, where upon a single toppling of a site the 
site looses two grains and gives one grain to each neighbor, except if the site is at the 
boundary where there is only one neighbor receiving a grain. See ^2] and na for more 
details on the abelian sandpile model in d = 1. 

The inhnite volume addition operator is dehned more precisely as follows: For 
r] E fl and x E with rj^x) = 1, we have axT] = f] + Cx (where ex{x) = 1 and ex{y) = 0 
otherwise), i.e., inactive becomes active, or, the height one at x simply changes to height 
two (and no other changes). For y E fl and x G Z with p(x) = 2 we look at the right 
- respectively at the left - of x to hnd the hrst site x+ (if that site does not exist we 
put x+ = oo), - respectively x~ (if that site does not exist we put x~ = —oo) - with 
= h(^^) = 1- We then dehne 



( 2 . 2 ) 


if both x+ and x^ exist. In words, upon adding one unit at x, the hrst sites at height 
one to the left (x“) and to the right of x (x+) become sites with height 2, and the site 
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which is the mirror image of x with respect to the middle of x^, and x~ becomes of 
height one, all other sites remain unaltered. We have to extend that definition to cases 
where one of the sites x^, x~ does not exist (i.e., when there is no site to the right or 
to the left of x having height one). That is done by taking the limit with “boundary 
condition 1”, i.e., if at least one of the x^ is infinite, then 


axViy) 


2 a x^ <y<x+ 
rjijj) otherwise 


(2.3) 


Remark that m is a special case of “addition” with “purely dissipative toppling”, 
i.e., upon toppling an active site two grains disappear (diagonal toppling matrix). In 
that sense combination of ax and 6x is the simplest example of combining two different 
toppling mechanisms (matrices) in one process. 


2.1 Construction 

Intuitively, our process is governed by two independent collections of Poisson processes, 
indexed by sites a; G Z, and independent for different sites. On the event 
times of we apply the addition operator, and on the event times of we “flip” 
the state, i.e., we apply 9x. We put the rate of the “sandpile-clocks” equal to a, and 
the rate of the “flip-clocks” equal to one. Formally, our process has as a generator on 
local functions / : O —> M, 

Lfiv) = - fiv))+ ~ 

:= aLsfiv) + Lf{v) (2.4) 

where Ls stands for “sandpile generator” and Lp for “flip-generator”. 

To show that there exists a Markov process with cadlag-paths corresponding to the 
Poisson process description above or to the formal generator Q , we use a monotonic¬ 
ity argument analogous to the one in M- We repeat the main steps and the minor 
modifications to be done here. First we define the action of a^, on as a “birth” if 
ri{x) = 1, and as an avalanche if r]{x) = 2, whereas the (identical) action of dx is (of 
course) also called a birth ii r][x) = 1, and a death if //(x) = 2. We can then split the 
formal generator in three parts: 

L = La -\- Li, -\- Ld (2-5) 

where 

Lafiv) = 

Lbfiv) = {l + a)J2x{v{x) = l){f{dxV) - fiv)) 

xGZ 

Ldfiv) = 5 ^x(^( 2 :) = 1)(/(M) -/(^)) (2.6) 

tcEZ 

Here, (x(-) denotes the indicator function). The construction is then as follows: 
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• Construct a process corresponding to La + Ld (only avalanches and deaths) on 
the set Qf oi conhgurations with a hnite number of sites with height 2. That is a 
(non-explosive) countable state space Markov chain on the set of hnite subsets of 
Z. Show by coupling that that process is monotone. The coupling is identical to 
that of 1121 for the avalanche events. For the deaths: we let two ones die together 
if possible, and otherwise independently. 

• Construct a process corresponding to La + Ld with births in a hnite interval, i.e., 
having generator 

n 

Lnfiv) = {La + Ld)f{ri) + xivix) = 1)(/(M) - fiv)) 

x=—n 

We construct that process once more as a countable state space Markov chain, 
and show that it is monotone. Its semigroup is denoted by Sn{t). Moreover, 
we have the following monotonicity as a function of the interval on which we allow 
births: for all f > 0, n G M, / monotone, p G 12/, 

{Ln{t)f){v) < {Sn+i{t)f){ri) 

• For general monotone / and p G 12 arbitrary: 

L{t)f{v) ■= sup sup Sn{t)f{r]) (2.7) 

nGN riGflf 


The process obtained by the above construction is called the SF-process (sand-hip 
process). We denote its path space measure starting from p by P,,. 

2.2 Basic properties 

Besides monotonicity, the SF-process has very similar “quasi-Feller” properties as the 
one-dimensional sandpile process of |I2|. In particular, we have the following analogue 
of Theorem 5.1 of ^21- Let us denote by 12' C 12 the conhgurations with an inhnite 
number of ones to the left and to the right of the origin. We then enumerate = 

{Xi{ri),i G Z} where Xoir}) := min{x > 0 : r]{x) = 1}, and the other Xj are in 
increasing order the sites where ri{x) = 1. The X, dehne the r^-dependent disjoint 
intervals Jj = (Xj_i(p), Xj(p)]. A function is called X-local if it depends on the heights 
rj{i) for i G Every local function is X-local, but a X-local function can be 

non-local, e.g. /(//) = jg bounded 1-local, but non-local. The idea is that 

the natural space to dehne the action of iterates of the generator is the set of X-local 
functions. That is made precise in the following dehnition and theorem. 

Definition 1 A configuration rj E fl' is called decent if 

1 

a{7]) = limsup— ^ |Wi(p) - Xi_i{ri) \ < oc (2.9) 

i=—n 

The set of decent configurations is denoted by Lldec- 
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Proposition 1 Let rj G fldec, f be bounded and N-local, then fort < 1/[4(1 + 0 ) 60 ( 77 )], 
the series converges absolutely and equals S(t)f{r]), where S(t) 

is the semigroup of the process defined above. In particular 


S{t)f{p) - /(77) 

t^o t 


Lfiv) 


i.e., L is the “pointwise generator” of the process. 


( 2 . 11 ) 


Proof. The same proof of |T2] can be used, if one notices that the extra “death” part of 
the generator can only split one of the intervals into smaller ones, by creating an extra 
1. This implies that if / is A^-local, then [/(oj?]) — /(t/)] = 0 for all 7 G Z \ 

Therefore Lf depends only on the heights in Iterating the argument, one 

sees that L^-f depends only on height in j, and one recovers the same estimate 

n /N+k \ 

ll(i")/|U < n E 2”(1 + “)"ll/llco (2.12) 

fc=0 Vi=0 / 

which gives the result of the theorem, by application of lemma 4.1 in ^21- • 

The following result, analogous to Corollary 6.1 in and to Proposition 3.1 in 
us. shows that the process is always non-Feller. 


Proposition 2 For all a > 0, the SF-process is non-Feller. 

Proof. We denote by 2 the maximal configuration rj = 2. We define the conhguration 

Vspec by 

. . f 1 if a: = 0 / , \ 

Vspecix) = W ,, . (2.14) 

I 2 otherwise 

Then one shows as in na that for fo{r]) = t]{0) 

= 2 (2.15) 

i.e., by the avalanche part of the dynamics, the isolated 1 is turned “immediately” into 
a 2. Therefore, the right limit of 17 * as t > 0, t > 0 is almost surely equal to 2 when we 
start from rjspec- 

This lack of right-continuity contradicts the Feller property. Indeed, if S{t) were a 
Feller-semigroup, then there would exist a uniformly dense set of continuous functions 
which are in the domain of the generator, i.e., for which 

lim = 11 ^/ 11 ^ (2.16) 

t^0,7>0 t 

However such a dense set of continuous functions contains a function / such that 

fiPspec) ^ /(2) (2.17) 

Combination of dnnD, ( 1 ^^ . (ITT7I) gives a contradiction. ■ 
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2.3 Stationary measure 

Denote by X the set of invariant probability measures of the SF-process dehned in the 
previous section, by S the set of translation invariant probability measures on D. By 
monotonicity of the process X is non-empty. In fact we have 

Theorem 1 For all a > 0, the SF-process is ergodic, i.e., X 
probability measures v on D, 

lim vSitd) = Pa 

t—5-00 

Moreover, 

• for a < 1, the density of sites with height one is given by 

j x(h(0) = l)dpa = (2.20) 

and Pa is a translation invariant measure which is mixing under translations, non¬ 
product and gives positive measure to all local events (i.e., has positive cylinders). 

• For a > 1, 

ha = ^2 (2.21) 

the Dirac measure concentrating on the maximal configuration rj = 2. Moreover 
for t > [log(Q; + 1) — log(Q; — l)]/2 and for every rj E fl, 

F.ipfiO) = 2 ) = 1 


= {pa} and for all initial 

(2.19) 


Proof. We start with the following lemma. 

Lemma 1 Let p be a probability measure on D that is mixing under spatial translations, 
with f x(h(0) = l)d/i = p > 0. Then we have 

a) Ift < p/[4(l-|-Q;)e], then pS(t) is mixing under spatial translations. In particular. 


lim [\S{t)[fT^g]-S{t)fS{t){Ta;g)\dp = 0 (2.23) 

fHoo J 

for all local functions f,g on D, where Tx denotes the spatial shift by x E Z. 
b) Lett\p) be the solution of 

(p+^)e-“'lP) + i^ = 0 (2,24) 

when it exists, otherwise put t'{p) = -|-cx). Then we have for all t < t\p). 


^ J viO)d{pS{t)) 

pit) ■= J X(h(0) = l)d{pS{t)) 


I Lp{0)d{pS{t)) (2.25) 

pe-2‘ + ^^(1 - (2.26) 
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Proof, a) Is exactly |T2], Lemma 6.1. 

b) Let to = mf{p/5(l + a)e,t'{p)}-, by a), ()2.25jl holds for t < to. Denote 

k~^{i,T]) = mi{j > 0-. r]{i + j) = 1} (2.27) 

k~{i,r]) = mi{j > 0 : r]{i — j) = 1} (2.28) 

We compute 

Lr]{0) = axiviO) = l)(fc+(l, r]) + k-{0,r]) + 1) + 3 - a - 2r]{0) (2.29) 

For u E S (see (6.77)), 

j X(h(0) = l)k-{0,r])du = j xivi^) = l)(^^(l,h) + = 1 

so that 

/ Lr,(0)d(^S(0) = « + 3 - 2 / ,md(,S(t)) = a - 1 + 2 /x(r,(0) = ms(t) (2.30) 
and hence 

^-M = -a + l-2m 

which gives for t < to. If to 7 ^ t'{p), then we can start the reasoning anew and 

iterate from to to p = min{p(to)/[5(l + a)e],t'(p(to))}, with new initial distribution 

pS(to). As a consequence, for 0 < t < to + p(t) is still given by (I2.26|l . etc. ■ 

By monotonicity, the limits 

Ui = lim 5iS{t) (2-31) 

t^OO 

for t = 1,2 exist and define invariant measures with ui < V 2 - Moreover, the process is 
totally ergodic (in the sense of HD], chapter 1, dehnition 1.9), i.e., X is a singleton if and 
only if Vi = V 2 - In that case vi = V 2 is also mixing under spatial translations, see pQ 
Theorem 1.4 ii) (indeed, the proof of the last theorem is a computation that does not 
involve the Feller property of the considered Markov semi-group, and relies on ()2.28ll L 
Let Xp denote the translation invariant product measure on D with Ap(77(0) = 1) = p. 
Then, using monotonicity and lemma ^ for t small enough, 

j X'())(0) = 1) d(S,S(t)) < j x(r,(0) = 1) d\S(t) = pe-” + 1^(1 - (2.32) 

Since p < 1 is arbitrary, we conclude 

j x(r,(0) = 1) d(S,S(,t)) < L^(l - e-“) (2.33) 

First consider a < 1. Starting from 5i, we have f'(p) = F(0) = oo, and we obtain 
from lemma H 

lim [ x(p(0) = 1) d{5iS{t)) = lim (^pe“2* + ^y-^(l - e’^*)^ = (2.34) 

t^OC J t^OO \ Z J Z 
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Similarly, for p close to one, t\p) = +cx), so that we can take the limit t > cx) in ()2.H2|1 . 
to obtain, using monotonicity once more. 


= lim [ xiviO) = l)d{5iS{t)) < lim [ x{v{0) = l)d{S 2 S{t)) < (2.35) 

Z t^OO J t^OO J Z 

Therefore, 

ni(r;(0) = l) = n2(r^(0) = l) = ^^ (2.36) 

and combining that with z/i < z/ 2 , gives vi = V 2 . 

To see that Pa is not a product measure, observe that because Pa ^ *5, it can be 
product only if p^ = Xp with p = (1 — a)/ 2 . Therefore, for all / local J LfdXp = 0. 
Indeed Xp concentrates in decent conhgurations, so we are allowed to use the generator 
by proposition ^ Notice that in the sandpile part of the dynamics two or more neigh¬ 
boring ones can never be created. Therefore, one computes the action of the sandpile 
part of the generator on the function Hn{p) = x(p(l) = ... = p{n) = 1) with n > 2, 
which gives after integration over the product measure 

/ C50 

LsH^dXp = -np^ - 2 *(1 - p)>-+i = -np^ - 2p^-\l - p) (2.37) 

i=0 

For the spin-flip part one has 

j LF{Hn)dXp = -np^ + np"-'(l - p) (2.38) 

Therefore, for the combined dynamics, the condition that Xp is stationary leads to 
j {aLs + LF){Hn)dXp = 0 = {—np + n{l — p) — nap — 2a(l — p)) (2.39) 


which gives (for n > 2 ) 


P 


n — 2 a 

2 n+ (n — 2 )a 


(2.40) 


Notice that equation is not valid for n = 1, because a single one can be created 

in the sandpile dynamics, so an extra term (1 — p) should be added for n = 1 . Since 
(CTHl should be valid for all n > 2, we obtain a contradiction. Hence, the invariant 
measure is indeed non-product. 

For a = 1, gives 


J x(h(0) = l)d{52S{t)) = 0 

Therefore 62 S(t) = 62 - On the other hand, ((H) gives 

lim [ x(p(0) = l)d{ 6 iS{t)) = 0 

t^ooj 
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Therefore 1^1 = 62 , hence we obtain ui = 1/2 = 62 - 
Finally, consider a > 1. Then we have 


Hence, for all £ > 0 there exists to < t^(l) such that for all t > to 

j x(?/(0) = l)d{5iS{t)) < e 

Therefore, z/i = limi^oo diS{t) = ^ 2 , and we conclude from the inequalities I'l = 62 > 1 ^ 2 ^ 
V 2 < ^2 that vi = V 2 = ^ 2 - 

Moreover, we obtain that from any initial distribution i/, the limiting measure 62 is 
reached in finite time <t'{l). I 

Remark 1 For a>l, we have f Lgd^a 7^ 0 for non constant g, since for the sandpile 
part f Lsgdjja = 0, whereas for the flip part f Lpgdp^a 0. Therefore, in that case the 
invariant measure cannot he found by solving J Lfdjj, = 0 for pi, which gives another 
argument for the non-Fellerian character of the SF-process. 

2.4 Generalization 

In this section we consider more general local perturbations of the sandpile generator. 
We show that the freezing phenomenon, i.e., having 82 as unique invariant measure for 
a large enough, and a non-trivial invariant measure for a small persists. 

More precisely, we consider a formal generator of the type 

L = OiLg T Lq (2.43) 

where Lq is the generator of a spinflip dynamics (i.e., with possibly configuration de¬ 
pendent rates): 


Lcfiv) = '^c{x,r]){f{9,,7]) - f{7])) (2.44) 

X 

where the flip-rates c{x, g) are supposed to be translation invariant, local and bounded 
from below. Therefore, 

m < c{x, g) < M (2-45) 

for some 0 < m < M < 00 independent of g. 

To define this process, we use a series expansion as in theorem [T] to define the 
semigroup. Remark that since we do not assume that Lq is the generator of a monotone 
process, the semigroup cannot be constructed by monotonicity. Instead, S{t)f{g) is 
defined by the series expansion as long as the configuration g is decent, and contrary 
to the monotone case, this cannot necessarily be extended to non-decent configurations 
such as the maximal configuration 2. 

We then have the following 
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Theorem 2 Consider the proeess with formal generator (EH- We have 


1. For a < m, 82 is not an invariant measure. In fact, if p. ^ S is an invariant 
measure, then 

777 — rv 

Fim = 1) > (2.47) 

2. If a > M, then for all p E S, with p{ri{0) = 1) > 0, pS(t) 82 as t ^ 00 . 
Therefore, 82 is the only possible invariant measure. 

Proof Let /i G 5 be such that p{r]{0) = 1) > 0. We denote pt = pS{t){ri{0) = 1). Then 
by the obvious generalization of lemma ^ we write 

^ = -a -2 j x{r]{^) = ^)c{ 0 ,r])d{pS{t)) + J c{ 0 ,r])d{pS{t)) (2.48) 

Therefore 

—a — 2Mpt + m < < —a — 2mpt + M (2.49) 

dt 

Hence, if a < m and pt < {m — a)/2M, 



so there can be no invariant measure p with p = p{p{If) = 1) < (m — a)/2M. That 
proves the hrst item of the theorem. For the second item, if a > M and if pt > 0, (j2.49j) 
gives 



and hence there cannot be an invariant measure with p = p{p{ 0 ) = 1 ) > 0 . ■ 


Remark 2 Even if pS{t) converges to 82 for all p E S with p{p{0) = 1) > 0, we cannot 
conclude that 82 is an invariant measure, because the process is not Feller. To see what 
can happen, consider the following example. Starting from a configuration p E Ft, we 
flip a 1 to 2 at rate 1 , independently for all lattice sites, and if the configuration is 2 , 
then we flip it at rate one to the minimal configuration p = 1, denoted by 1. This is 
a non-Fellerian process with pS(t) 82 for all translation invariant measures p 7 ^ 82 , 
but clearly, 82 is not invariant. In fact, the process has no invariant measure. 

If Lc is the generator of a monotone process, then more precise results can be obtained. 
For that case we will stick to an explicit example where once more an explicit closed 
equation for the density can be obtained. More precisely, we consider the flip rates 

c{x,p) = l-7^(h)(/x-i(h) + /x+i(h)) (2.51) 

where 

fxiv) = 1 - 27 ( 7 ( 3 ;) = 1 ) 
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These rates correspond to the standard Glauber choice for 

7 = ^ tanh(2/3) e 

where (3 denotes the inverse temperature (here without any meaning except for an 
effective coupling constant). We then have the following analogue of theorem ^ 

Theorem 3 For the process with formal generator and rates (irni) we have 

1. For a < Uc = I — 27 , there exists a unique non-trivial invariant measure pa with 

A«('((0) = 1) = - ^) (2-53) 

2. For a > ac , S 2 is the unique invariant measure. 

Proof. Since the rates (EUl) satisfy definition 2.1 of chapter 3 in ma, the process with 
generator Lq is monotone. Therefore, we can construct the process with generator 
f|2.43|l by monotonicity as in section 12. IL where we replace the coupling for the birth 
and death part by basic coupling. In particular the thus obtained generalized SF-process 
is monotone. For /i G 5 with /i( 7 ( 0 ) = 1) > 0 we obtain 

^ = -« + (1-2a)(1 - 27) (2.54) 

and from that equation, combined with monotonicity we can proceed as in the proof of 
theorem ^ ■ 

Remark 3 As one would expect intuitively, the critical value etc is decreasing in 7 , i.e., 
the freezing is enhanced by stronger coupling. 

Another simple choice in which we explicitly see the effect on ac is obtained by adding 
a bias to the spin flip. Then, becomes 

c{x,r]) = 1 - Kf^{ri) = (1 - K)xiv{x) = 2) + (1 + K)x{vix) = 1) 

and a similar calculation yields the same result but with a critical value that now equals 

ac = 1 — K. 

3 Adding “anti-additions” to the sandpile process 

3.1 The anti-sandpile model 

In words, the anti-sandpile process is a process where grains are removed from a con¬ 
figuration 7 G G, and afterwards (if necessary) the configuration is stabilized instanta¬ 
neously by reversed topplings. 


1 1 
2 ’ 2 
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We first define the finite-volume process. If after removing grains the height is 
zero at one or more sites x G [—iV, iV], then the conhguration stabilizes by a sequence 
of reversed topplings. Upon a reversed toppling of a site x G [—N, iV] the site gains 
two grains and each of its neighbors (in [—iV, iV]) looses one grain. This means that 
in a reversed toppling, the boundary sites act as a source (instead of a sink in the 
ordinary toppling rule). The anti-addition operator aj, is then dehned as the stable 
result of the subtraction of one unit at site x and performing reversed topplings until 
the conhguration is stable (i.e. height everywhere 1 or 2) again. 

Remark 4 The anti-addition operator should not he confused with the inverse of the 
addition operator. In fact, if rj is recurrent, and air] is recurrent, then alrj = af^i^rj), 
but alirj) need not he recurrent if rj is. 

In hnite volume, the generator of the anti-sandpile process is given by 

N 

= E (“i - -f) P'2) 

x=-N 

where / denotes the identity operator. Remark that aj, = Oa^O and 

= QLQ (3.3) 

where 0 is “global spinhip”, i.e., 

0/(h) = f{0v) (3.4) 

with 6 {r]){x) = {r]{x) 1) mod 2. Therefore the extension of the process generated by 

Li to inhnite volume is immediate. Its semigroup is given by 

S{t)^ = QS{t)e (3.5) 

where S(t) is the semigroup of the sandpile process. 

In inhnite volume, the “anti-addition operator” aj, is then dehned via 

alrj = 6ax0{r]) 


Similarly, we introduce 


l^{x, r]) = k^{x, Or]) 


and the intervals Ji{r]) = liiOrj). 


3.2 The SA process 

We now dehne the SA-process (i.e. “sandpile -|- anti-sandpile”) as the process associated 
to the formal generator 

Lap = o:L -t- /3L^ (3.6) 

where L = Ls = - I)- 
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This process is constructed as follows: we define the semigroup acting on local 
functions via the series expansion of proposition ^ This gives the hnite dimensional 
distributions, and hence dehnes a unique Markov process starting from decent conhgu- 
rations, where decent means here that both rj and Orj are decent in the sense of dehnition 
m We call the thus dehned process the “SA-process”. 

We then have the following. 

Proposition 3 The SA-process is monotone. As a consequence, it can be defined start¬ 
ing from any initial configuration. 

Proof. In ini we constructed a generator for a coupling of the sandpile process which 
preserves the order. The idea of this coupling is simply that if rj < ^, then for each site 
j having height two in ^ which by an addition at some site z G Z could be turned into 
a one, in rj either the height of j is one or there exists a unique site x{j, rj, having 
height two in rj such that addition at that site creates in // a site of height one. 

Let us call the (formal) generator of this coupling. Remark now that if ?7 < 
then of course Q{ri) > Q{^), and for all / monotone, 0(/) is also monotone. 

Therefore, the coupling with generator 

{LJ){p,O = c^iLyMO + (^iL's{Qf)m,0O (3.8) 

dehnes a coupling that preserves the order. 

This proves monotonicity. The consequence is clear since every conhguration can 
be written as an increasing (or decreasing) limit of decent conhgurations. ■ 


3.3 Stationary measures for the SA process 

We then have the following analogue of theorem ^ 

Theorem 4 LetT be the set of invariant measures for the process with generator (ESI). 
Then we have 

• For a < j3 

I={ 6 ,} 

• For a > (3, 

T = {^ 2 } 

• For a = 

{51,62} 

Proof. We compute as before, starting from a translation invariant measure /i concen¬ 
trating on decent conhgurations: 

j LsxiviO) = f)dT = -1 
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and hence since x(h(0) = 1) + x(h(0) = 2) = 1: 

j = +1 

Hence, starting from an initial measure /i on which is translation invariant, mixing 
and concentrates on decent conhgurations, we obtain, using once more the notation 
Pt = I Sit){xiv{0) = l))d/i: 

^ = - paL + /3L')(XO(0) = l))<i(fiS(()) = (/3 - a) (3.10) 

Of course this equation is only valid as long as 0 < (/3 — a)t < 1. It expresses that the 
density of ones simply decreases or increases linearly until no ones are present, resp. all 
sites are of height one. 

Notice that we used here the analogue of proposition d which in this case implies 
that one can use the generator as in the Feller case as long as acting on local functions 
and integrated over measures that have a non-zero density of sites having height two 
and of sites having height one. 

Starting from this equation, one concludes that for a > j3 (and similarly for a < (3), 
there can be no other invariant measure (which is also translation invariant) than 62 
(resp. (5i). We then deduce the hrst two statements of the theorem along the same lines 
as in theorem ^ using monotonicity ( by proposition IHl). 

For the last statement, use that if a = (3, 


dpt 

dt 


= 0 


Therefore, using standard arguments based on monotonicity, we see that and 62 are 
invariant measures. ■ 


Remark 5 For the case a = (3 we have 



i.e., the density is a conserved quantity. An open question here is whether in that case 
for each density there exists a stationary (in time) and ergodic (under translations) 
measure with that density or whether the only extremal invariant measures are {(5i, 62 }- 

In that case, we can however say the following: 

1. If p is translation invariant, invariant for the dynamics, and with density 0 < 
p < 1, then /i© is also invariant for the dynamics. Indeed, for any decent 
function f, 0/ is also decent, and J{Ls + QLsQ)fdp = 0 is equivalent to 
J{Ls + eLse){Qf)d{pQ) = 0 . 
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2. The product measure \p with density Xp{ri{t)) = T) = p is not invariant. Indeed, 
we can proceed as in the proof of theorem{Jl and compute, for Hniji) = x(^(l) = 
... = r]{n) = 1) with n >2, 



(3.12) 
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